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Abstract 

We use lattice techniques to calculate the continuum string tensions of SU (N) gauge theories in 
2 + 1 dimensions. We attempt to control all systematic errors at a level that allows us to perform 
a precise test of the analytic prediction of Karabali, Kim and Nair. We find that their prediction 
is within 3% of our values for all N and that the discrepancy decreases with increasing A^. When 
we extrapolate our results to = cxd we find that there remains a discrepancy of ~ 1%, which is a 
convincing ~ 6a effect. Thus, while the Karabali-Nair analysis is remarkably accurate at = oo, 
it is not exact. 

PACS numbers: 
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I. INTRODUCTION 



Finding an analytic solution to confinement and the mass spectrum of four dimensional 
SU{N) Yang-Mills gauge theories continues to be a challenge for theoretical physics. It is 
in the large- iV limit [ij that the relation between gauge field theories and string theory is 
most natural and this relation is strengthened by the AdS/CFT correspondence. Indeed 
a solution to the problem of confinement may shed light on string theories with moderately 
strong coupling. 

In lower dimensions, gauge theories are of interest for the same theoretical reasons. In D = 
1 + 1 all dynamical degrees of freedom can be removed by gauge fixing, and the theory can 
be analytically solved, with linear confinement arising trivially from the Coulomb potential. 
Moving up one dimension to D = 2 -|- 1, makes the theory much more complicated. Gluons 
become dynamical, a linearly confining potential appears to be dynamically generated, and 
the theory appears to be as analytically intractable as in the four dimensional case. This 
is unfortunate since the D = 2 + 1 theory is also phenomenologically interesting: through 
dimensional reduction, it approximates the high temperature limit of the four- dimensional 
theory. Indeed, an analytic solution in D = 2 + 1 dimensions would be a significant step 
forward, and would perhaps bring us closer to a solution of the D = 3 + 1 case. 

While an exact solution to D = 2 + 1 gauge theories is not available, Karabali and 
Nair have claimed to obtain a very good approximation through studying the continuum 
Hamiltonian, and expressing it in terms of color-singlet fields [3]. Truncating the Schrodinger 
equation for their ground state functional, they obtained the following prediction for the 
tension a of the string that binds distant static sources in the fundamental representation: 



Remarkably, this prediction turns out to be within ~ 3% of the lattice calculated values 
0, Q for all values of N. 

If one replaces these sources with ones in a general representation, TZ, then the analysis 
of j3] predicts that the tension a-ji of the string between these sources, obeys 



where C-ji is the quadratic Casimir of the representation TZ. In Eq. (11.21) there are no 
restrictions on TZ, and so this clearly cannot be exact in general: in D = 2 + 1 gluons 
screen color charge and will, for example, break the string if TZ has trivial A^-ality. However 
in the limit iV — > oo screening vanishes and, moreover, the lattice calculations jsl, 0] show 
a decreasing discrepancy with Eq. (11.11) as N increases. Indeed the extrapolated N = oo 
lattice value is only about ~ 1% below the prediction of Eq. (11.11) . At this level of accuracy 
there are a number of systematic errors that were inadequately controlled in these lattice 
calculations, and so this leaves open the tantalising possibility that the calculation of [3] 
might become exact at N = oo. The purpose of the present paper is to provide a lattice 
calculation in which all the systematic errors are controlled at a level that allows us to test 
this possibility. 

One of these systematic errors has to do with corrections to the string energy as a function 
of length. This is an interesting theoretical subject of its own, having to do with the nature 
of the effective string theory that describes the confining flux tube. We expand in detail on 




(1.1) 



(1.2) 
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this issue in a companion publication [7]] and here will only quote a few relevant results. A 
related issue is how the string tension varies as a function of the representation of the flux. 
We extend a similar control of systematic errors to these calculations in a second companion 
paper Q where we test how well the Casimir scaling in Eq. (11.21) is satisfied - in particular 
by stable fc-strings. 

In the next section we describe the basic lattice setup. We then describe the methodology 
used in the calculation of the string tension, what are the important systematic errors, and 
how we control them. We then provide our results, extrapolate them to the continuum limit, 
and then extrapolate to = cxd. We finish with our conclusions. 

Some of our preliminary results, both from this paper and from Q , have been presented 




II. LATTICE SETUP 

The theory is defined on a discretized periodic Euclidean three dimensional space-time 
lattice, with spacing a and, typically, with L^Lt sites. The Euclidean path integral is given 
by 

Z = J DUexp{-pSw), (2.1) 

where f3 is the dimensionless lattice coupling, and is related to the dimensionful coupling g"^ 
by 

2N 

lim/?=— . (2.2) 

In the large- A^ limit, the 't Hooft coupling A = g'^N is kept fixed, and so we must scale 
(3 = 2N'^/X cx A^^ in order to keep the lattice spacing fixed (up to 0(1/A^^) corrections). 
The action we choose to use is the standard Wilson action 



Sw — X! 



1 - —ReTi Up 
N 



(2.3) 



where P is a lattice plaquette index, and Up is the plaquette variable obtained by multiplying 
link variables along the circumference of a fundamental plaquette. We calculate observables 
by performing Monte-Carlo simulation of Eq. (12. ip . in which we use a mixture of Kennedy- 
Pendelton heat bath and over-relaxation steps for all the SU{2) subgroups of SU{N). 



III. METHODOLOGY 

To obtain the string tension, a, we calculate the energy of the lightest flux tube that 
winds around one of the spatial tori. Extracting the mass from the correlation function is 
the first area in which we need to control systematic errors, as described below. 

In the confining phase such a winding flux tube cannot break. Since there are no sources 
here (in contrast to Wilson loop calculations of the static potential) there are no extraneous 
contributions to the energy (such as the Coulomb potential) and one can hope to obtain a 
string description for all lengths / of the flux tube. There is of course a smallest possible 
length /c = l/T'c, where Tc is the deconfining temperature, below which there are no winding 



flux tubes. However for A^ > 4 the transition is first order [lOl] and so for larger A^ we can hope 
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to have a string description for any / > Ic- Such a description should become particularly 
simple at = oo where mixing and decay vanish. We have performed a careful study of the 
way the string energy depends on its length, which will be published elsewhere Here we 
use those results to bound the theoretical uncertainties in extracting the asymptotic string 
tension from the string energy, so controlling the associated systematic errors. 

There are also systematic errors in extrapolating to the continuum a = limit and, 
subsequently, to the N = oo limit. These will be discussed below as well. 



A. Extracting string masses from correlation functions 

We calculate the energy of the winding flux tube from the correlator of suitably smeared 
p = Polyakov loops that wind around a spatial torus. This is a standard technique [sl, 0] 
with the smearing/blocking designed to enhance the projection of our operators onto the 
ground states. (We use a scheme that is the obvious dimensional reduction of the one in 
We calculate with several blocking levels and construct the full correlation matrix. From 
this we obtain best estimates for the ground and excited string states using a variational 
method applied to the transfer matrix T = e~°'^ - again a standard technique jl, 0, |Tl|, 12 



In practice, our best operator for the string ground state has an overlap ~ 99% onto that 
state so that the normalised 'ground state' correlation function satisfies 

= (1- |e|)exp{-Mo(Ot} + |ei|exp{-Mi(/)t} + ... ; ^ [e^l = |e| ~ 0.01 (3.1) 

i 

where Mq, Mi are the ground and first excited state string energies. (Since our time-torus is 
finite, we use cosh fits rather than simple exponentials, although in practice we use Lt large 
enough for any contributions around the 'back' of the torus to be negligible.) To extract 
Mq from this correlator one can fit with a single exponential for t > Iq, discarding the lower 
t values so that a statistically acceptable fit is obtained. This is a reasonable approach 
and one followed in [l, 0]. However it neglects the systematic error arising from the fact 
that there is certainly some excited state contribution as demonstrated, for example, by the 
fact that one cannot obtain a good fit with a single exponential from t = 0. To control 
this systematic error we also perform fits with two exponentials, with a fixed mass M* for 
the excited state, resulting in a mass Mo(M*) for the ground state. Typically Mo(M*) is 
smallest when M* is as small as possible, i.e. M* = Mi, and is largest when M* = oo, i.e. 
effectively a single-exponential fit. So the true value typically satisfies: 

Mo (Ml) < M^'""" < Mo(oo). (3.2) 

From here on, we refer to the single-cosh fitting procedure by 'S', and to the double-cosh 
fitting procedure, by 'D', and add these as superscripts to any relevant results (such as 
masses, string tensions, etc.). Consequently we shall have two continuum string tensions , 
and a^, that bracket the true string tension 

< ^truc < (3 3) 
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B. Extracting string tensions from string energies 



From the ground state string energy, Mq{1), we need to extract the tension cr. Taking 
into account the Liischer term [1^], a is given by 

Moil) vr , , 

^ = ^ + Z77 + 77 • (3-4) 



/ 6/2 ' 

In practice using -\/cr/ > 3 one can expect the neglected 0(1//"^) corrections to be smalL 
However they represent another systematic error that needs to be controlled. To do this we 
shall use the results of our study in 0], where we have calculated Mq{1) as a function of / 
for N = 3,4, 6, 8. We have done so in the range 1.3 — 1.6 < ly/a < 3 — 6.2 and for different 
lattice spacings. We find that our results can be well encompassed by 

(Moil)y TT 0.2(1) 

V al ) 3{^iy i^lY' ^ ■ ' 

Here the first two terms on the right hand side constitute the exact Nambu-Goto string 
prediction. It is believed that when one expands Mo{l)/al in inverse powers of y/al the 
first two corrections are universal 13, l3| and equal to those for the Nambu-Goto string. 



(This is consistent with our numerical calculations in [3].) The next possible correction term 
corresponds to the last term on the right of Eq. (13.51) . We observe that the fitted coefficient 
is much smaller than the 0(1) coefficient characteristic of the other terms. This shows that 
Nambu-Goto provides a remarkably accurate description of the ground state winding string 
energy for all possible lengths. 

The calculations of Mq{1) in this paper are performed for l^/a > 3. We extract the 
corresponding values of cr using Eq. (13.51) . For such / the contribution of the correction 
terms that are additional to the Liischer correction is in fact almost negligible, as was 
assumed in earher calculations jl, 0]. However now we are able to control the accuracy of 
that assumption. 



C. Extrapolation to the continuum limit 



To extrapolate to the continuum limit we need to choose a theoretically motivated fitting 
ansatz for the way ay/ a depends on the bare lattice coupling f3. From Eq. (12. 2p it is clear 
that 



P 



lim 



(3.6) 



The leading perturbative correction to this relation will be 0{l/f3). (Note that we do not 
have here a ratio of physical quantities for which the leading correction would be 0(a2).) 
Lattice perturbation theory is notoriously ill-behaved and to reduce the higher order correc- 



tions we replace (3 by the mean-field improved coupling [15 



/?MF = X (TlUp) 



(3.7) 




the continuum extrapolation was performed with a leading 0(1//3mf) correction. 
The values of ay/ a from the coarser values of a typically did not lie on the fit, and were 



5 



excluded. Although this is a sensible procedure, it ignores the (small) systematic error due 
to the neglect of higher order corrections in 1//?mf- Here we will control this error by fitting 
with an additional correction 



By comparing these fits with linear fits where we constrain 02 = 0, we shall have an estimate 
of the effect of higher order corrections. Clearly such a strategy is only possible where the 
calculations are of sufficient range and accuracy, as in the present paper. 

It is important to note that this procedure is not without its ambiguities. The expansion 
in Eq. (13. 8p is a weak coupling one which is functionally incorrect in the strong coupling 
region. If our fit includes one or more points in the strong coupling region, it is these 
points that may well determine our estimate of the coefficient 02 in Eq. (13.81) . in which 
case the estimate will be unrelated to its actual value. This problem is exacerbated by the 
fact that unlike the case in D = 3 + 1 (for > 5), the separation between strong and 
weak coupling does not involve a clear-cut first order transition, but rather a cross-over 
that turns into a smooth phase transition at = cxd very much like the Gross- Witten 
transition in D = 1 -(- 1 [17|. This crossover peak increases with A^ and lies in the range 
(3/2N'^ ~ 0.40 — 0.45 [l6|, and so in our fits we shall avoid using any values obtained on the 
strong coupling side of this peak. 

D. Extrapolation to = cxd 

The continuum value of y/a / g'^N is expected to have a finite limit, with leading corrections 
of 0{1/N'^). Hence linear fits in 1/A^^ can be used to extrapolate to A^ = 00, as in 0]. 
To control the neglected higher order corrections to the linear fit we shall also perform fits 
using the more general form 



and compare the results to those of linear fits (62 = 0). There is no reason to believe that 
this expansion becomes functionally incorrect at small A^, and we shall use it all the way 
down to A^ = 2. 

IV. RESULTS 

In Fig. [1] we plot the values of ay/aPuF/i'^'N'^), as obtained from both S and D fits, 
versus N^/ Pmf for our SU(4) and SU(6) calculations. We display fits that are quadratic, as 
in Eq. (13. 8p . The 0(l//5^) correction is always positive, except possibly for SU(2), so the 
quadratic fit leads to a significantly higher value than the linear one in the continuum limit. 
Despite this fact, we see from the figure that the extrapolated string tensions are still lower 
than the values predicted by Karabali, Kim, and Nair in Eq. (II. ip . 

We give the results of our continuum fits in Tables HHITl where we also list for comparison 
the Karabali-Kim-Nair predictions. 




(3.8) 




(3.9) 
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0.2 r-r 




FIG. 1: The dimensionless quantity (3mfj^ — ^ as a function of the improved inverse couphng 
couphng 1//3mf for = 4, 6. The error bars at 1/Pmf = denote the result of the continuum 
extrapolation, while the horizontal bars denote the values predicted by Karabali, Kim, and Nair 

a. 





iV = 2 


= 3 


iV = 4 


iV = 5 


iV = 6 


iV = 8 


af 


-0.0133(7) 


-0.0158(4) 


-0.0186(9) 


-0.0175(9) 


-0.0204(8) 


-0.0186(8) 


a? 


-0.0122(12) 


-0.0152(9) 


-0.0168(19) 


-0.0160(14) 


-0.0200(11) 


-0.0186(10) 


af 


-0.0004(3) 


+0.0005(2) 


+0.0023(6) 


+0.0018(6) 


+0.0040(5) 


+0.0025(5) 


"2 > 


-0.0009(5) 


+0.0004(3) 


+0.0012(12) 


+0.0010(8) 


+0.0039(7) 


+0.0025(6) 



TABLE I: The parameters ai_2 in the fit Eq. (13. 8p . which are obtained for the string tensions a^'^. 
The superscripts S, D denote the way we fit the correlation function, and bracket the actual string 
tension (see Section UlI Ap . 



Since we are interested in the accuracy of the Karabai-Kim-Nair (KKN) prediction with 
increasing A^, we define the ratio r of that prediction to our lattice values 

We now fit with the form 

= + ^ + If. ("-2) 
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Type of fit 


Source of data 


N = 2 


iV = 3 


iV = 4 


iV = 5 


iV = 6 


iV = 8 


Linear 




0.1678(1) 


0.1837(2) 


0.1892(1) 


0.1917(2) 


0.1932(2) 


0.1948(2) 


Linear 




0.1675(3) 


0.1831(3) 


0.1886(3) 


0.1910(4) 


0.1927(3) 


0.1943(3) 


Quadratic 




0.1675(3) 


0.1839(2) 


0.1902(3) 


0.1924(3) 


0.1944(3) 


0.1955(3) 


Quadratic 




0.1666(6) 


0.1832(4) 


0.1893(7) 


0.1915(5) 


0.1939(4) 


0.1951(4) 


KKN prediction 


0.17275 


0.18806 


0.19314 


0.19544 


0.19668 


0.19791 



TABLE II: Our continuum values of ^fa jg^N and the predictions of Karabali, Kim, and Nair 
(KKN) H- We present results from linear and quadratic extrapolations to the continuum [0,2 = 
or 02 as a free fit parameter). 

and present the results for linear fits (with = 0) and quadratic fits (with C2 as a free 
parameter) in Table UTTl below. (The number of degrees of freedom for all fits was d.o.f = 3.) 
Finally, we show in Fig.[2]the ratio r plotted against with S and D linear fits. As for the 
continuum extrapolation, the effect of the higher order correction is to lift the N = 00 value 
towards the Karabali-Nair value, albeit not quite far enough to achieve a perfect match. 



Type of fit 


Source of data 


^00 


Cl 


C2 


xVd.o.f. 


Linear 




0.9902(12) 


-0.215(32) 




0.84 


Linear 




0.9878(17) 


-0.245(56) 




0.64 


Quadratic 




0.9908(14) 


-0.268(58) 


0.41(22) 


0.98 


Quadratic 




0.9886(21) 


-0.317(103) 


0.52(40) 


0.57 



TABLE III: The extrapolation of r to the large- limit. 



V. SUMMARY 

We have calculated the tensions of strings in the fundamental representation of SU{N) 
gauge theories in 2+1 dimensions. Our immediate goal was to test the prediction of Karabali- 
Nair in Eq. (11. ip . particularly at N = 00, where the screening effects that are clearly not 
incorporated in that scheme, vanish. Since earlier lattice calculations had already shown 
that any discrepancy was no more than a few percent, our calculation needed to control a 
number of previously neglected systematic errors that are small but could be significant at 
the ~ 1% level. 

In this paper we described how we controlled the following errors. Firstly the contribution 
of excited states to our variationally selected ground state correlators, from which we extract 
the energy of the ground state winding flux loop. Secondly higher order string corrections in 
the relationship between this ground state energy and the asymptotic string tension. (Using 
the results of our companion publication [3] on the effective string theory describing winding 
flux loops.) Thirdly, higher order corrections in the continuum extrapolation, and fourthly, 
higher order corrections in the extrapolation in to = 00. 

Our final results are similar to the ones of the older work which assumed that the 
systematic errors that we control here, are negligible. We find that this assumption is, as it 
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- data 


O 


1 i noQ r fit 
— Ml lcd.1 ML 




- data 




- linear fit 




1/N^ 

FIG. 2: The ratio r between the prediction of Eq. (II. ip and our data, as a function of 1/N'^. The 
error bar at = denotes the hnear extrapolation to the N = oo hmit. 



happens, essentially correct. Our string tensions are 3% — 1% smaller than the prediction 
of Eq. (11. ip . and a discrepancy persists at = oo, where our result is 



= 0.1975 ±0.0002 - 0.0005. (5.1) 

\^ / Lattice 

Here the first error is statistical, and the second comes from the difference 6(7 = as — ctd > 0. 
This error can only lower the string tension, away from Eq. (11.11) . Consequently, our result 
is lower by 0.98% - 1.2% than Eq. ffTTTj) 



KKN 



0.199471 . . . 



(5.2) 



This difference while small is statistically significant at a compelling 8 — 5.4 sigma level 
(depending on the details of the fit). 

While it is clear that the leading term in the scheme of j3] is not exact at N = oo, our 
results show that it is astonishingly accurate. 
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